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Abstract 

The problem of a compact binary system whose components move on circular orbits is addressed 
using two different approximation techniques in general relativity. The post-Newtonian (PN) 
approximation involves an expansion in powers of v/c <C 1, and is most appropriate for small orbital 
velocities v. The perturbative self-force (SF) analysis requires an extreme mass ratio m\jmi <C 1 
for the components of the binary A particular coordinate-invariant observable is determined as 
a function of the orbital frequency of the system using these two different approximations. The 
post-Newtonian calculation is pushed up to the third post-Newtonian (3PN) order. It involves 
the metric generated by two point particles and evaluated at the location of one of the particles. 
We regularize the divergent self-field of the particle by means of dimensional regularization. We 
show that the poles oc (d — 3)" 1 appearing in dimensional regularization at the 3PN order cancel 
out from the final gauge invariant observable. The 3PN analytical result, through first order in 
the mass ratio, and the numerical SF calculation are found to agree well. The consistency of this 
cross cultural comparison confirms the soundness of both approximations in describing compact 
binary systems. In particular, it provides an independent test of the very different regularization 
procedures invoked in the two approximation schemes. 
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I. INTRODUCTION 



A. Motivation 

The detection and analysis of the gravitational radiation from black hole binaries by 
the ground-based LIGO- VIRGO and space-based LISA observatories requires very accurate 
theoretical predictions, for use as gravitational wave templates (H. There are two main 
approximation schemes available for performing such calculations in general relativity: (i) 
The post-Newtonian expansion, well suited to describe the inspiralling phase of arbitrary 
mass ratio compact binaries in the slow motion and weak field regime (c _1 = v/c <C l), 1 
and (ii) the perturbation-based self-force approach, which gives an accurate description of 
extreme mass ratio binaries (q = m\jm2 <C 1) even in the strong field regime. 

For the moment the post-Newtonian (PN) templates for compact binary inspiral have 
been developed to 3.5PN order in the phase jl-HJ and 3PN order in the amplitude [f| @] (see 
j8[ for a review). These are suitable for the inspiral of two neutron stars in the frequency 
bandwidth of LIGO and VIRGO detectors. For detection of black hole binaries (with higher 
masses) the PN templates have to be matched with full numerical simulations for the merger 
phase and the ringdown of the final black hole. The matching between the PN approximation 
and numerical relativity has turned out to be very successful [9|, 1 101) . 

On the other hand, gravitational self-force (SF) analysis |lll-ll5| is expected to provide 
templates for extreme mass ratio inspirals (EMRIs) anticipated to be present in the LISA 
frequency bandwidth. SF analysis is a natural extension of first order perturbation theory, 



and the latter has a long history of comparisons with post-Newtonian analysis [16H23J. SF 
analysis, itself, is just now mature enough to present some limited comparisons with PN 
analysis, but it is not yet ready for template generation. 

In this paper we shall compare the PN and SF analyses in their common domain of 
validity, that of the slow motion weak field regime of an extreme mass ratio bi nary (see 
illustration of various methods in Fig. [1]). The problem was tackled by Detweiler |24(, who 
computed numerically within the SF a certain gauge invariant quantity, defined by (11. 3p 
below for an extreme mass ratio binary, and compared it with the 2PN prediction extracted 



from existing PN results [25[. Here we shall go one step further, and extend the comparison 
up to 3PN order. This will require an improvement in the numerical resolution of the SF 
calculation in order to distinguish more accurately the 3PN self-force from the self-force at 
higher PN orders. However, our primary difficulty is that the PN results for the metric have 
not previously been available at 3PN order, and will have to be carefully derived. We shall 
demonstrate an excellent agreement between the extreme mass ratio case C 1) of the 
analytical 3PN result and the numerical SF result. 



B. Method 

Let us consider a system of two (non-spinning) compact objects with masses m\ and rri2, 
and moving on slowly inspiralling quasi-circular orbits. In the PN analysis, let ni\ and m>2 



By a slight abuse of notation we denote by c _1 the standard PN estimate, where c is the speed of light. 
As usual we refer to nPN as the order equivalent to terms 0(c~ 2n ) in the equations of motion beyond the 
Newtonian acceleration. 
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FIG. 1: Different analytical approximation schemes and numerical techniques are used to study black hole 
binaries, depending on the mass ratio mi/7712 and the orbital velocity v 2 ~ Gm/rvi-, where m = mi + ni2- 
The post-Newtonian theory and black hole perturbation theory can be compared in the slow motion regime 
(v <C c equivalent to T\i 3> Gm/c 2 for circular orbits) of an extreme mass ratio (mi -C 1112) binary. 



be arbitrary; in the SF analysis, further assume that mi <C mi- We can then call m\ the 
"particle" , and the "black hole" . 

Self-force analysis shows that the dissipative parts of the self-force for a circular orbit 
are the t and tp components. These result in a loss of energy and angular momentum from 
the small mass at the same precise rate as energy and angular momentum are radiated 



awav|24| . In addition, earlier perturbative calculations of energy and angular momentum 



fluxes [I6H23J for this situation show them to be equivalent to the results of the PN analysis 
in their common domain of validity. Hence, by invoking an argument of energy and angular 
momentum balance, we know that the PN results also agree with the dissipative parts of 
the SF in their domain of common validity, and further comparison can reveal nothing new. 

For our PN-SF comparison, we shall thus neglect the dissipative, radiation-reaction force 
responsible for the inspiral, and restrict ourselves to the conservative part of the dynamics. 
In PN theory this means neglecting the dissipative radiation-reaction force at 2.5PN and 
3.5PN orders, and considering only the conservative dynamics at the even-parity 1PN, 2PN 
and 3PN orders. This clean separation between conservative even-parity and dissipative 
odd-parity PN terms is correct up to 3.5PN order. 2 In SF theory there is also a clean split 
between the dissipative and conservative parts of the self-force. This split is particularly 
transparent for a quasi-circular orbit, where the r component is the only non- vanishing 



2 However, this split merges at 4PN order, since at that approximation arises a contribution of the radiation- 
reaction force, which originates from gravitational wave tails propagating to infinity [26]. 
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component of the conservative self-force. 

Henceforth, the orbits of both masses are assumed to be and to remain circular, because 
we are ignoring the dissipative radiation-reaction effects. For our comparison we require two 
physical quantities which are precisely defined in the context of each of our approximation 
schemes. The orbital frequency Q of the circular orbit as measured by a distant observer is 
one such quantity. The second requires further explanation. 

With circular orbits and no dissipation, the geometry has a helical Killing vector field k a . 
A Killing vector is only defined up to an overall constant factor. In our case k a extends out 
to a large distance where the geometry is essentially flat. There Jc^<9 Q = dt + d v in any 



natural coordinate system which respects the helical symmetry [27|. We let this equality 
define the overall constant factor, thereby specifying the Killing vector field uniquely. 

An observer moving with the particle m 1 , while orbiting the black hole 7712, would detect 
no change in the local geometry. Thus the four-velocity uf of the particle is tangent to the 
Killing vector k a evaluated at the location of the particle, which we denote by kf. A second 
physical quantity is then defined as the constant of proportionality, call it uj, between these 
two vectors, namely 

u a x =u\kl. (1.1) 

The four- velocity of the particle is normalized so that (p Q!i a) 1 tt"ztf = — 1; {g a p) 1 is the regu- 
larized metric at the particle's location, whereas the metric itself is formally singular at the 
particle mj in both PN and SF approaches. The gauge invariant quantity uj is thus given 
by: 



' -(a^*?*?) 1 = (-(a*)!*?*?) -1 ' 2 . (1.2) 



It is important to note that this quantity is precisely defined in both PN and SF frameworks, 
and it does not depend upon the choice of coordinates or upon the choice of perturbative 
gauge; however, it very definitely depends upon using a valid method of regularization. 
Furthermore, for any coordinate system uj has a pleasant physical interpretation as being 
the rate of change of time at a large distance, with respect to the proper time on the particle 
mi, and it could in principle be measured by a redshift experiment as described in |24j . 

If we happen to choose a convenient coordinate system where k a d a = dt+fl everywhere, 
then in particular k\ = 1, and thus uj = u\, the t component of the four velocity of 
jri]. The Killing vector on the particle is then kf = uf/u\, and simply reduces to the 
particle's ordinary post-Newtonian coordinate velocity v®/c. In such a coordinate system, 
the description of the invariant quantity we are thus considering is 



a 8\ -1/2 
I t~ \ v f v i 



< =A= ( -GMx-^rJ • (1.3) 

In the PN calculation we shall evaluate uj using a particular harmonic coordinate system. 
We shall make no restriction on the mass ratio q = m^mj, but shall eventually compute 
the small mass ratio limit q <C 1 for comparison with the SF result. 

The regularized metric (g a p)i is defined with very different prescriptions in the SF and 
PN approaches. Both analyses require subtle treatment of singular fields at the location 
of the masses. Subtracting away the infinite part of a field while carefully preserving the 
part which is desired is always a delicate task. Our comparison will rely on the principle 
of the physical equivalence of the regularized SF and PN metrics, at least in the vicinity 
of the particle, i.e. that they are isometric — they differ by a coordinate transformation. 
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In fact the cross cultural comparison of the invariant uj is a test of the isometry of the 
two regularized metrics and is, thus, a test of the two independent (and very different) 
regularization procedures in use. 

In the SF prescription, the regularized metric reads 



9ap(x) = 9a P (x) + Kp{x) , (1.4) 

where g a p denotes the background Schwarzschild metric of the black hole, and where the 
"Regular" perturbation is smooth in a neighborhood of the particle, and follows from the 
Detweiler- Whiting prescription [l3| for removing the infinite part of the field, as described 
below in Sec. Ill Bl In particular the metric fll.4p is regular at the particle's position y", and 
we simply have 

affix = flSG/i)- (1.5) 

In the perturbative SF analysis we are only working through first order in q = 777,1/7712, and at 
that level of approximation = 0(q). Then uj can be computed accurately to the same 



perturbative order and compares well with the post-Newtonian result to 2PN order [24 
The regularized 2PN metric is known [25| , and therefore the comparison is straightforward. 

In the present paper we shall obtain the 3PN regularized metric which will be the core 
of our calculation, and will be partly based on existing computations of the equations of 



motion at 3PN order using Hadamard [28| and dimensional [29J regularizations. Using an 



iterative PN procedure, one first considers the post-Newtonian metric #fg(x, t) at any field 
point outside the particle, in a coordinate system x a = {ct, x 1 }. That metric is generated by 
the two particles, and includes both regular and singular contributions around each particle. 
Then we compute the PN regularized metric at the location of the particle by taking the 
limit when x — > yi(t), where yi(t) is the particle's trajectory. In 3 spatial dimensions, that 
limit is singular. In order to treat the infinite part of the field, we extend the computation 
in d spatial dimensions, following the prescription of dimensional regularization, which is 
based on an analytic continuation (AC) in the dimension d viewed as a complex number. 
Considering the analytic continuation in a neighborhood of e = d — 3— > 0, we define 

(^^ACflim^x,*)!. (1.6) 



x-tyi 



The limit e — > does not exist in general due to the presence of poles oc e' 1 occurring at 
3PN order; we compute the singular Laurent expansion when e — > 0, and we shall see that 
the poles disappear from the final gauge invariant results. Previous work on equations of 
motion and radiation field of compact binaries has shown that dimensional regularization is 
a powerful regularization method in a PN context. In particular this regularization is free 
of the a mbig uities plaguing the Hadamard regularization at the third post-Newtonian order 
I, & 0, M ■ 



The plan of this paper is as follows: Sec. [TT] is devoted to an overview of the SF formalism. 
The circular geodesies of the perturbed Schwarzschild geometry are described in Sec. Ill Al 
where we also give an explicitly gauge invariant relationship between Q and uj for the particle 
mi. We use the mode-sum regularization procedure of Barack and Ori [U 32| to perform 



the delicate subtraction of the singular field h^a from the retarded metric perturbation h r ^. 
We give a brief description of our application of this process in Sec. Ill Bl In Sec. Ill CI we 
describe some of the details of the numerical analysis which yields our value for as a 
function of f2, and provide a brief discussion of the numerical determination of the 3PN 
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effect on uf . Most of the details concerning the 3PN calculation are presented in Sec. IIHI 
We focus mainly on the issues regarding our implementation of the dimensional self-field 
regularization which is described in Sees. IIII Al and IIII B[ The post-Newtonian results are 
presented in Sec. [TV] We give the fully-hedged regularized 3PN metric in Sec. IIV A[ and 
present our final result for uj in Sec. IIVBI We finally investigate the small mass ratio limit 
q <C 1 of the post-Newtonian result, and compare with the self-force calculation in Sec. [Vj 
Two appendices provide further details on the PN calculation: An alternative derivation 
using the Hadamard regularization is discussed in Appendix El and the choice of the center- 
of-mass frame and the reduction to quasi-circular orbits at 3PN order within dimensional 
regularization are investivated in Appendix [B] 



II. SELF-FORCE OVERVIEW 

Previously we described the truly coordinate and perturbative-gauge independent prop- 
erties of Q and uj . In this section we use Schwarzschild coordinates, and we refer to 
"gauge invariance" as a property which holds within the restricted class of gauges for which 
k a d a = dt+Q d v is a helical Killing vector. In all other respects, the gauge choice is arbitrary. 
With this assumption, no generality is lost, and a great deal of simplicity is gained. 

The regularized metric perturbation h^a = h^g—h^g is the difference between the retarded 
metric perturbation h™g and the singular field h^g- A Hadamard expansion of Green's 



functions in curved spacetime provides an expansion for h^g 13j |. In a neighborhood of the 
particle with a special, locally-inertial coordinate system, h^g appears as the mi/r part 3 of 
the particle's Schwarzschild metric along with its tidal distortion caused by the background 



geometry of the large black hole. Details of the expansion are given in Sec. 6.1 of (33j. The 



special locally inertial coordinates for a circular geodesic in the Schwarzschild metric are 



given as functions of the Schwarzschild coordinates in Appendix B of [34 



A. Circular geodesies of the perturbed Schwarzschild geometry 



The effect of the gravitational self-force is most easily described as having mi move along 
a geodesic of the regularized metric g aB + h^ g . We are interested in circular orbits and let 
u a be the four- velocity of mi. 4 This differs from the four- velocity u a of a geodesic of the 
straight Schwarzschild geometry at the same radial coordinate r by an amount of 0(q). 
Recall that we are describing perturbation analysis with q <C 1, therefore h^g = 0(q), and 
all equations in this section necessarily hold only through first order in q. 

It is straightforward to determine the components of the geodesic equation for the metric 



9aB 



24 



and then to find the components of the four-velocity u a of mi when it is 
in a circular orbit at Schwarzschild radius r. We reiterate that the four-velocity is to be 
normalized with respect to g a g + h^g rather than g a8 , and that h^g is assumed to respect 



the symmetry of the helical Killing vector. In this case we have 



t\2 



3mo 



1 + vTvPh*:* 



-u a u 



8a uR 

o r n aB 



(2.1a) 



3 In all of Sec. U we set G = c = 1. 

4 Since we are clearly interested in the motion of the small particle mi, we remove the index 1 from u". 
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r(r — 3m 2 ) 



m 2 {l + u a vPh%) 1 

7 \ — u u p d r h„ R 

r(r-2m 2 ) 2 aP 



(2.1b) 



A consequence of these relations is that the orbital frequency of mi in a circular orbit about 
a perturbed Schwarzschild black hole of mass m 2 is, through first order in the perturbation, 
given by 

The angular frequency Q is a physical observable and is independent of the gauge choice. 
However the perturbed Schwarzschild metric does not have spherical symmetry, and the 
radius of the orbit r is not an observable and does depend upon the gauge choice. That 
is to say, an infinitesimal coordinate transformation of 0(q) might change u a vfid r h^^. But 
if it does, then it will also change the radius r of the orbit in just such a way that Q 2 as 
determined from (12. 2p remains unchanged. Both u l = u T and = Qu T are gauge invariant 
as well. 

Our principle interest is in the relationship between Q and u T , which we now establish 
directly using (12. lap and (I2.2p . writing all equations through first order. First, we can 
rearrange (12.21) to get: 

^' = (m2 Q f + (!^) 2 (! - ^i) gwrtW*,) . (2.3) 

Next, we take the cube root of both sides and expand on the right-hand-side (RHS) to 
obtain: 

The second term on the RHS of (12.41) is already first order in q. Thus, in the first two 
bracketed expressions in this second term, we can replace m 2 /r by the leading approximation 
to m 2 /r from just the first term on the RHS of (12 .4p . giving: 



^ = (m 2 fi) 2 / 3 + Vl - 3(m 2 Q)^ (L u a u%h^\ . (2.5) 



Following [24| , we next introduce the gauge invariant measure of the orbital radius 
Now we use this in its second form and substitute back into (12.51) : 

Multiplying overall by —3 and adding 1 to both sides before dividing through, we find: 

' 1 - T - u a u p d r hl^] = o \ /p . (2.8) 



3m 2 /r V 2 r a(i J 1 - 3m 2 /R n ' 
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This is exactly what we need in ( 12. lap in order to establish a first order, gauge invariant, al- 
gebraic relationship between u T (to which u l evaluates in our gauge) and Rq (or equivalently 
Q), namely: 

W 2 =(l-^y\l + u°tfh%). (2.9) 

The lowest order term on the RHS is identical to what is obtained for a circular geodesic of 
the unperturbed Schwarzschild metric. Indeed, recall that the Schwarzschild part of u T is 
known exactly as ^s chw = (1 — 3m 2 / Rq)' 1 ^ 2 . Thus, if we write 

u T = ul hw + qul F + 0(q 2 ), (2.10) 

the first order term in (I2.9P gives: 

9U&= K 1- ^r) V ^*' (2 - n) 

which is 0(q), and contains the effect of the "gravitational self-force" on the relationship 
between u T and Q, even though it bears little resemblance to a force. We shall henceforth 
focus our attention on the calculation of the combination u a u^h?:a. 



B. Mode sum regularization 

Both the retarded metric perturbation h™p and the singular field h^g are singular at vn\. 
However, we actually determine h™p by using the inherent symmetries of the problem to 
separate variables and to decompose the components of h™g in terms of tensor spherical 

harmonics. Each £, m component h 1 ^ is then finite and determined using a standard 
numerical differential equation solver. Only the sum over modes diverges. 

For our problem, we treat the divergence of the singular field h^g in a related manner. 
The singular behavior is represented in the known expansion of h^g about the particle, and 
is also amenable to a decomposition in terms of spherical harmonics. This procedure is 
stylistically quite similar to the expansion of the Coulomb field of a point charge, displaced 
from the origin, in terms of spherical harmonics centered on the origin; this results in the 
coefficients being proportional to either l/r e+l or r , depending upon whether the field point 
is inside or outside the charge. In SF analysis, the spherical harmonic coefficients determine 
the regularization parameters of h^g- 



Following the original prescription of Barack and Ori [31|, |35| and extending it as in (34 



we first perform the sum over m for the retarded field at the particle 

e 

fh r *f ] = u a u^h x ^ m) . (2.12) 



u vr 

ULIJ 

m=-l 



Then we use the recognition that the decomposition of the singular field is of the form 
u u p h\ = B + - — — - + — — — — — — + 



£+1/2 (2i-l)(2e + 3) (2£-3)(2£-l)(2£ + 3)(2£ + 5) 
+ 0(r 6 ), (2.13) 
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where B, C, D, E\ (and the subsequent E2, E3, etc) are regular izat ion parameters. The 
particular £ dependence of the coefficients accompanying the parameters D and E n is related 
to the expansion of (1 - cos#) n+1 / 2 in terms of Legendre polynomials P^(cos#); details are 
derived and described in Appendix D of (34j . 

The regular field at the particle is finally given by 

^% = E(^C W - « afi/, > ( 2 - 14 ) 

and the sum is guaranteed to be convergent as long as B and C are known. In practice, the 
regularization parameters are difficult to determine. For our problem it is known analytically 
that C = and 



3m 2 



2m<2 



where 2 -Pi is a hypergeometric function, and r is the Schwarzschild radial coordinate of 
the circular orbit. This knowledge of B and C, but not D, implies that the sum in ( 12. 14ft 
converges as 1 /£. To increase the rate of convergence, we augment our knowledge of B and C 
by numerically determining further regularization parameters (IH : We use the fact that the 
behavior of u a vPfi^ , evaluated at the particle, must match u a u^h^ as given in (I2.13P for 
large i. This allows us to fit the numerical data to determine the additional regularization 
parameters D and E n up to, say, E%. Knowledge of these additional parameters results in 
a sum which converges as l/£ 9 . In our numerical work we typically fit for three or four 
extra parameters. We calculate up to £ = 40, fit in the range £ = 13 — 40, and then sum to 
i — > +00, with errors at the full level of our calculational precision. 



In earlier work [24j, |27|, |36| the accuracy of the numerical integration used was adequate 
for the purposes then at hand. For the comparisons presented here it became obvious 
that we should investigate pushing our integration procedure to enable us to obtain the 
highest precision practicable. By adjusting the effective step size as £ changed, 5 we found 
that it was possible to achieve this without encountering any other numerical difficulties 
(such as an unreasonable accumulation of round-off error). Subsequent fitting, to obtain 
the numerical determination of the higher order regularization parameters D, • ■ ■ , E3 as 
described above, allowed us to reduce residuals to the level of the computational precision 
which had controlled our integration procedure. Monte Carlo calculations based on these 
residuals gave us systematic estimates of the errors to associate with our fit parameters. 
Using these, we find relative errors of order 10 -13 in u^ F (the loss in precision being due to 
the regularization). The corresponding results are presented in Table HI 

After the regularization procedure is complete we have in hand u a u l3 h^ j3 , and hence Wg F , 
for an orbit at a given radius R^. At this point, we have solved our self-force problem - 
we have found the effect, «g F , of the self-force on u T for a specific Q. 



5 We used an integration procedure that contained an adjustable parameter, e, which controlled the precision 
of the numerical result. We chose e small enough so that further reduction would not cause relative changes 
in the result larger than 1CP 15 . 
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Rn/m 2 




-a-BuR 1 
u a u p h2p/q 




T 
n SF 


200 


-0 


0100252390238679 


-0 


00505064245513028 


220 


-0 


00911174844278219 


-0 


00458725834137915 


240 


-0 


00835083080996084 


-0 


00420175898117037 


260 


-0 


00770720725494635 


-0 


00387603022007156 


280 


-0 


00715569723937482 


-0 


00359717107497568 


300 


-0 


00667784659538770 


-0 


00335574417643231 


320 


-0 


00625982212277844 


-0 


00314468649077390 


340 


-0 


00589105041112645 


-0 


00295860680303681 


360 


-0 


00556331104384481 


-0 


00279331869895365 


380 


-0 


00527011654983391 


-0 


00264552181684313 


400 


-0 


00500627861027562 


-0 


00251257921031088 


A90 


n 

— u 


uu^ i o i oyoooooyooz 


n 




440 


-0 


00455064124356486 


-0 


00228311728867935 


460 


-0 


00435257068802445 


-0 


00218341682793388 


480 


-0 


00417102337921533 


-0 


00209205962311231 


500 


-0 


00400401451882955 


-0 


00200804044413982 



TABLE I: Summary of the gravitational self- force effects for a variety of radii Rq. Approximately 
13 digits are believed to be accurate. 



C. Post-Newtonian fit of «g F 

The improved quality of the data in Tabled fed directly into the next stage, that of fitting 
iigp as a function of Q (or Rq) to determine the higher order post-Newtonian coefficients. 
In order to proceed to our post-Newtonian fit of Mg F , we introduce a special notation for the 
convenient gauge invariant PN parameter defined in (12.61) . which is 0(c~ 2 ) and reads 

y = (m 2 Qf 3 = ^ . (2.16) 
The post-Newtonian expansion of the self-force effect given in (12. lip was determined explic- 



itly up to 2PN order in [24J , and found to be 

i& = -y - 2y 2 - 5y 3 + Cf£ N y 4 + 0(y 5 ) , (2.17) 



where C|p N represents the 3PN coefficient (unknown at the time of 24|). and higher-order 
4PN terms are neglected. The author of Ref. {24j also performed a numerical fit of iig F to 
the polynomial ( I2.17P in order to determine the numerical value of C|p N . He expected that 
the post-Newtonian derivation of this coefficient would be unavailable for some time, and 
the numerical fit was done in a cursory fashion using a range in Rq from 20 m 2 to 50 rri2, not 
generally optimal for PN comparison. It is now time to improve upon that early analysis. 

The process of fitting terms in the £-sum for the regularization parameters is relatively 
easy. Convergence in the £-sum increases by two orders with each additional regularization 
parameter, and is very rapid. By contrast, the low order PN series for Ug F is effectively a 
power series in 1/Rq and is relatively slowly convergent. Moreover, if we tried to fit higher 



10 



terms, we could rapidly encounter the situation where, for some large Rq, contributions 
would be below our error estimates, while for smaller Rq, the same contributions would still 
be significant. This situation complicated both the choice of the range of Rq over which we 
could effectively fit, and the choice of the number of additional PN coefficients we should 
use to improve the characterization of our available data, consistent with the error estimates 
we had previously established. For this paper, we settled on a compromise, namely, we used 
values of ifo generally in the range 200 m 2 to 500 m 2 and, surprisingly, somewhere between 
3 and 6 additional PN coefficients. Following these procedures, our numerical determination 
of the 3PN coefficient in the self-force effect upon u T was found to be 

C 3 S | N = -27.677 ± 0.005 . (2.18) 

If we do not use a sufficient number of additional coefficients, our fitting procedure will 
compensate by systematically attempting to approximate the missing terms by adjusting 
the fitting parameters we do use. This effect, which determines the error term in (12.181) . 
turns out to dominate the random error from our numerically determined data points. We 
can estimate this effect by our choice of the range of data and the number of coefficients used 
in the fitting process. This proves to be adequate for our comparison while also suggesting 



that further work is warranted. In particular, we shall show in separate work |37J that the 



PN expansion involves in higher orders some logarithmic terms, and that the prior knowledge 



of the coefficients of the logarithms appearing at 4PN and 5PN orders (computed in [37 
from PN theory) will improve very much the accuracy of the PN fit to the SF result. 6 This 
study is beyond our current scope, but will be extensively reported in j^ - " 



III. POST-NEWTONIAN CALCULATION 



In this section, our aim is to compute the 3PN regularized metric (jl.6p by direct post- 
Newtonian iteration of the Einstein field equations in the case of singular point mass sources. 
Previous work on the 3PN equations of motion and radiation field of point particles |H, 0, H^, 



30| has shown that the appropriate regular izat ion to remove the infinite self-field of point 



particles in this context is dimensional regularization |38|, [39 . 

In the dimensional regularization (DR) scheme, we look for the solution of the Einstein 
field equations in D = d + 1 space-time dimensions, with a matter source made of point par- 
ticles. We treat the space dimension as an arbitrary complex number, rfeC, and interpret 
any intermediate formula in the PN iteration of those equations by analytic continuation in 
d. Then we analytically continue d down to the value of interest (namely 3), posing 

d = 3 + e. (3.1) 

In most of the calculations we neglect terms of order e or higher, i.e. we retain the finite 
part and the eventual poles. 

Defining the gravitational field variable h a/3 = yf—g g al3 — t] al3 , 7 and adopting the har- 
monic coordinate condition d^h^ = 0, we can write the "relaxed" Einstein field equations 



6 Accordingly, the 0(y 5 ) symbol for remainders in Eq. (|2.17[) and similar equations below should rather be 
understood as the Landau o(y 4 ) symbol. 

7 Here g a ^ is the contravariant metric, inverse of the covariant metric g a p of determinant g = det(c/ Q ^), 
and r] a " = diag(— 1, 1, 1, 1) represents an auxiliary Minkowski metric in Cartesian coordinates. 
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in the form of ordinary d'Alembert equations, namely 



where □ = rf' v d iJ d v is the /?a£-spacetime d'Alembertian operator in D space-time dimen- 
sions. The gravitational source term A af5 in ( 13. 2 p is a functional of h py and its first and 
second space-time derivatives, and reads as 

- g^guadph^d^ - g fi »g va d p h aa d t Ji'"> + g^g pa d p h^d a h^ 

+ i(2<7<V v - <7 a V) (g P *ge. - ■^- [ g* e g P x)d li h' m d v h m . (3.3) 

Note the explicit dependence on the space dimension d of this expression. The matter 
stress-energy tensor T a/3 will be composed of Dirac delta-functions in d dimensions, say 
§( d ) [x _ y a ] ; where x is the field point and y a are the source points labeled by a. Finally the 
d- dimensional gravitational constant is related to the usual Newton constant G by 

G {d) =Gf , (3.4) 

where £q denotes the characteristic length associated with dimensional regularization. We 
shall check in Sec. IIVI that this length scale never appears in the final 3-dimensional result. 



A. Post-Newtonian metric in d dimensions 



The 3PN metric is given in expanded form for general matter sources in terms of some 
"elementary" retarded potentials (sometimes called near- zone potentials) V, V$, K, Wy, B4, 
X, Zij, Yi and T, which were introduced in Ref. 28| for 3 dimensions [see Eqs. (3.24) there] 
and generalized to d dimensions in Ref. 29j. All these potentials have a finite non-zero 
post-Newtonian limit when c — > +00 and parameterize the successive PN approximations. 
Although this decomposition in terms of near- zone potentials is convenient, such potentials 
have no physical meaning by themselves. Let us first define the combination 



r , 4 A 16T 
K+ — + — 



(3.5) 



Then the 3PN metric components can be written in the rather compact form 2^ 



So P o N 



_ -2V/c 2 



SViVi 32RiVi 



CP 



0{c 



10\ 



(3.6a) 



PN 



_g (d-2)c 2 



I ( d ~ 1 V 

+ 2 \ J^2c^ 



+ 



8Ri 



16 

77 



+ 0(c 



(3.6b) 



This particular exponentiated form is to be consistently recxpandcd at 3PN order. 
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9 



PN 



g (d-2)c 2 



16 



W + 



2(d-2) 



+ 0(c- 



(3.6c) 



The successive PN truncations of the field equations f)3.2p — (13.3H give us the equations satisfied 
by all the above potentials up to 3PN order. We conveniently define from the components of 
the matter stress-energy tensor T a/3 the following density, current density, and stress density 



a 

era 



2 (d - 2)T 00 + T 



d-1 

rpOt 

! 

c 



where T n = 5ijT % K The leading-order potentials in the metric obey 
UV = - 4nG {d) a , 

DK = — AnG {d) aV , 



r &kk 



\(d-l 



2\d-2 



diVdiV. 



(3.7a) 

(3.7b) 
(3.7c) 

(3.8a) 
(3.8b) 
(3.8c) 

(3.8d) 



These potentials evidently include many PN corrections. The potentials V and Vi have a 
compact support (i.e. their source is localized on the isolated matter system) and will admit 
a finite limit when e — > without any pole. With the exception of the potential K which 
has also a compact support, 9 all other potentials have, in addition to a compact-support 
part, a non-compact support contribution, such as that generated by the term oc diVdjV 
in the source of Wij. This is the non-compact support piece which is the most delicate to 
compute because it typically generates some poles oc \/e at the 3PN order. The d'Alembert 
equations satisfied by all higher-order PN potentials, whose sources are made of non-linear 
combinations of lower-order potentials, are reported here for completeness: 



URi = - 



AirGW f5-d Tr d-1 

-d^2\rr VcTi -^2 

d-1 „ „„„ d(d-l) 



nx 



d-1 

47rG {d) 



dnVdiVk 
Vo, 



4(d-2) 



dtVdiV 



(3.9a) 



d-2 



+ 2 



d-3 



d-1 



1 (d-1 



+ W ij d ij V + 2V i d t d i V^ 



d-3 
d-2 



VdfV 



al T + K 



d{d 



A(d-2) 2 



{d t Vf - 2diVj djVi 



(3.9b) 



9 Actually the compact-support potential K does not contribute to the present calculation. Indeed, it will 
always be multiplied by a factor e = d — 3, and being compact docs not generate any pole; so it does not 
exist in 3 dimensions. 
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uz 



'J 



□T 



AuG {d) 
d-2 



V an - 5, 



ij 



d-2 



d 



d-2 



dtViidfiV + diVkdjV* 



d k Vi d k Vj - 2d k V(i dj)V k 
d(d-l 



Si ; 



8(d-2y 5 - {dtVy 



d k V m (d k V m — d m V k ) 



d-2 
(d-l)(d-3) 



1 



2{d-2) 



d-2 
& kk Vi 



2(d-2) 2 
(5-d)(d 



(3.9c) 



4(d-2) 2 



aVVi 



1 * 1 

2 CkW ik + -o-jfcVfc 



d- 3 
(d^ 



(J; 



v- 



5 -d 



A' 



1 (d 



d-2 
d(d 



d k v + <9 fc v; - d k w tl d t v k 



T=^ dkV9ikk 4(d-2) 2 



1 V k d t Vd k V- i {d 7 ^ VdtVdiV 



1 



+ 



d-2 
(d- l)(d-3) 

47rG (d) 



&AaK + 



d-2 
d(d- l)(d- 3) 
4(d-2) 3 



Kd-2) ; 
7fi^ + 2V fc d k dtVi 



(dtVdiK + diVd t K) 



3 



1 

12 



d- 3 



d-3 



d-2 



(3.9d) 
aX 



aV A - 



1 /d-3 



2 V d 



aVK + 



(5-d)(d-3) 
2(d- l)(d-2) 



OiViV 



d-3 f- 
d — 1 



d-3 
2(d- 2) 2 
1 /d 



er,;A 



4(d-2) 2 



d-2 

2 



d 



d-2 



VidjVidjV 



8 Vd-2 



^ 2 ^V + 



d(d 



Kd-2) ; 



v(d t vy 



1 (dV? ^-^-V vd 2 R d(d-l)(d-3) 
2 (<TO " 4(d-2) 2 4(d-2)3 dtVdtK 

(d-l)(d-3) 2 ^-3 l/d-3 
4(d-2) 2 J ^ ^ ~ d^2 Vi mK ~ 2 \d—2 



(3.9e) 



Many of the latter potentials have already been computed for compact binary systems, 
and we shall extensively use these results from 28|, |29|. Notably, all the compact-support 
potentials such as V and Vi, and all the compact-support parts of other potentials, have been 
computed for any field point x, and then at the source point yi following the regularization. 
However, the most difficult non-compact support potentials such as X and T could not be 
computed at any field point x, and were regularized directly on the particle's world-line. 
Since for the equations of motion we needed only the gradients of these potentials, only the 
gradients were regularized on the particle, yielding the results for (0jX)(yi) and (<9jT)(yi) 
needed in the equations of motion. However the 3PN metric requires the values of the 
potentials themselves regularized on the particles, i.e. X(yi) and T(yi). For the present 
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work we have therefore to compute, using the tools developed in |28|, |29| , the difficult non- 
linear potentials -X"(yi) and T(yi), and especially the non-compact support parts therein. 
Unfortunately, the potential X is always the most tricky to compute, because its source 
involves the cubically-non-linear and non-compact-support term Wij dijV, and it has to be 
evaluated at relative 1PN order. 

In this calculation we also meet a new difficulty with respect to the computation of the 
3PN equations of motion. Indeed, we find that the potential X is divergent because of 
the bound of the Poisson-like integral at infinity. 10 Thus the potential X develops an IR 
divergence, in addition to the UV divergence due to the singular nature of the source and 
which is cured by dimensional regularization. The IR divergence is a particular case of 
the well-known divergence of Poisson integrals in the PN expansion for general (regular) 
sources, linked to the fact that the PN expansion is a singular perturbation expansion, with 
coefficients typically blowing up at spatial infinity. The IR divergence will be treated in 
Sec. IIIIB2I " 

The 3PN metric 13.61 is valid for a general isolated matter system, and we apply it to the 
case of a system of N point-particles with "Schwarzschild" masses m a and without spins 
(here a = 1 , ■ • • , N) . In this case we have 

cr(x, t) = Y,K 5 {d) [x - y a (t)} , (3.10a) 
<7i(x, t) = $> a < 6& [x - y a (t)} , (3.10b) 
^•(x, t) = A*« < < S {d) [x - y.(t)] , (3.10c) 

a 

where 6^ denotes the Dirac density in d spatial dimensions such that J d d x5^(x) = 1. We 
defined the effective time-varying masses of the particles by 

» a (t) = ma = , (3.11) 

y{ggai3)(ya,t)v% v^/c 2 

together with p, a = ^ [d - 2 + v^/c 2 ] /i a . 



B. Dimensional regularization of Poisson integrals 

In the PN approximation we break the hyperbolic d'Alembertian operator □ in Eqs. (13. 8p - 
(13. 9p into the elliptic Laplacian operator A and the small PN retardation term c~ 2 df, which 
is put in the RHS of the equation and iterated. Neglecting the radiation-reaction effects, this 
means that we solve the d'Alembert equations by means of the symmetric Green function 

□sVm = A" 1 + ^A" 2 «9 t 2 + 0(c" 4 ) . (3.12) 

We consider only the 1PN retardation because the potentials X and T which are the only 
ones to be computed are to be evaluated at 1PN order at most. We are thus led to define 
the dimensional regularization of Poisson or twice-iterated Poisson integrals. 



However the other potential T, which is merely Newtonian, is convergent at infinity. 
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Let F(x) be the generic form of the functions representing the PN potentials in d di- 
mensions. For simplicity we shall treat only the case of the non-compact support terms. 
Compact support potentials or compact part of potentials, such as V or the first term in the 



source of do not generate poles in d dimensions and were dealt with in Ref. [29j. Also, 
we consider only Poisson integrals and refer to [29| for the procedure for iterated Poisson 
integrals. So we want to compute a typical Poisson potential 

P(x') = A" 1 [F(x)] = - A / d<x . (3.13) 



Ait J |x — x\ 

We employ the Green function u = k |x| 2 ~ d of the Laplace operator in (^-dimensions, satis- 
fying An = — Att8^(x.), where 

k = -A^l (3.14) 

is defined from the standard Eulerian gamma function. 11 Furthermore we want to evaluate 
the Poisson integral -P(x') on one of the singular points, say x' = y a . As we shall see the 
Poisson potential we have to deal with will not only be divergent on the singularities but 
also at infinity, i.e. when the source point r = |x| — > +oo. To delineate these problems we 
introduce a constant radius 1Z > 0, and split the Poisson potential into a near-zone integral 
corresponding to r < TZ, and a far-zone one such that r > TZ: 

P(x') = P<(x') + P>(x% (3.15) 

The near-zone integral P< will contain the local or ultra-violet (UV) singularities, due to 
the singular nature of the point-mass source, while the far-zone integral P> will have an 
infra-red (IR) divergence at infinity, which is actually a general feature of the PN expansion 
for any post-Newtonian source. 



1. UV divergence 

The UV divergences will be dealt with using dimensional regularization. Non-compact 
support terms are generated by a generic function P(x) which extends to all space. For all 
needed cases we can assume that F is smooth everywhere except at the singular points y a , 
around which it admits a singular expansion in powers of r a = |x — y tt | when r a — > 0, of the 
type 

P qi 

f w = EE r « +q£ f™(*°) + °W > ( 3 - 16 ) 

p=po q=qo 

for any P G N. The coefficients a fp,q(^a) are functions of the unit direction n a = (x — y a )/r a , 
and depend on the dimension through e = d — 3, and also on the DR length scale £q. 12 The 
powers of r a are of the type p + qe, where p and q are relative integers (p, q G Z) with values 
limited as indicated. The singular expansion (I3.16P will yield some UV-type divergence of 
the Poisson potential ( 13.1 3p . Relying on analytic continuation, we can evaluate F at the 



11 The constant k tends to 1 when d — > 3, and was formerly denoted k in Ref. [29(. 

12 More precisely, the coefficients depend on £ as a f P , q oc l^ qi ' , as can be seen from the expansion (|3.16|) . 
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location of particle a simply by taking the limit x — > y a . Indeed, we can check that the 
dimension d can always be chosen such that F is non-singular in this limit. Thus, 

F(y ) = Ac[lim F(x)l , (3.17) 

and we may then consider the expansion when e — > 0. From now on the analytic continuation 
process will be implicitly assumed without indication. 

The near-zone part of the Poisson integral of the generic function F outside the singu- 
larities is defined by 

in which the upper bound of the integral is set at the intermediate radius TZ. The singular 
behavior of this integral at the source points, i.e. when x = y tt , is automatically taken care 
of by dimensional continuation down to d — 3. Next we evaluate the integral at the singular 
point x' = y a itself. The result is easy in DR, 13 as we are allowed to simply replace x' by 
y a into (I3.18p . Thus, 

k f ^J(x) 
'r<n 

which is the main result of DR, as applied to UV divergences. 

In practical calculations we are interested in the three-dimensional limit, so we perform 
the Laurent expansion of the previous result when e = d — 3 — > 0. As we know from previous 
work [29|, the expression (I3.19P is finite for any of the non-compact potentials up to 2.5PN 
order, but will develop a simple pole oc 1/e at the 3PN order. The poles correspond to the 
occurrence of logarithmic divergences in the three-dimensional calculation (28| . and are in 
fact associated with our particular coordinate choice. Indeed, similar calculations performed 



P < W = -5L fc ^' (319) 



at the 3PN level in ADM-like coordinates within DR are pole-free [30[. The information 
we shall need is the pole part followed by the associated finite part when e — > 0; we shall 
usually leave aside the remainder term 0(e). lA We thus consider the expansion 

P<(y a ) = - £ Pt\y a ) + Pf(y«) + 0{e) , (3.20) 

and we look for the pole-part P< ^(y a ) and finite-part P<\y a ) coefficients. Beware that 
our terminology is slightly misleading, because we shall conveniently include in the pole part 
P< l \ya) some dependence in e, which will of course be thought as being expanded when 
e —7- up to first order in e, therefore yielding a finite contribution O(e ) to be added to 
the finite part P< (y ). Combining previous results in Sec. IV of [29J, we find that the pole 
part is explicitly given by 



1 A 



e 

Q=Q0 



q a q+ b^a e=o l - Vr «fe / b 



(3.21) 



13 This is in contrast with the difficult formulation necessary in Hadamard's rcgularization; see Refs. 281.129 1. 

14 Although those remainders 0(e) present in Newtonian terms will be kept because they might get multiplied 
by some poles e -1 at the 3PN order, therefore yielding finite contributions at 3PN order. 
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The first term is the contribution of the singularity a which is clearly singled out, while the 
second term comes from all the other singularities b ^ a. The bracket notation in (13.2 ip 
refers to the angular average performed in d dimensions, i.e. 

(U - / %^ fM , (3-22) 
where the volume Qd-i of the (d — l)-dimensional sphere is given by 

a ^ =m = w=rv (3 ' 23) 

We observe that the pole part (I3.2ip depends on the behavior of the function F only for x 
in a neighborhood of the singularities, through the singular expansion coefficients Q / M with 
p ^ —2. The result for the pole part heavily relies on DR, and depends on the DR scale 
£q introduced in (I3.4p . However, as it is "localized" on the singularities, the pole part is 
independent of the radius 1Z. 

On the contrary, the finite part coefficient P< (y ) depends on all the "bulk" of the 
integration outside the particle's world-lines and in particular on the radius 1Z. This finite 
part essentially corresponds to what we would naively compute in 3 dimensions, i.e. by 
simply replacing d = 3 into (13.1 9ft ; the result, however, would be ill-defined as it stands. In 
DR the finite part coefficient in (I3.2(jp has a precise meaning, and we find that it agrees with 



the so-called Hadamard partie finie integral |40|, |41 



Pf(y a ) = -^Pf, / d 3 x^, (3.24) 



An 



r<Tl 



where F^> is the function F computed with e — 0. Here "Pf^ " stands for the partie finie 
which depends on the arbitrary scale £q playing here the role of the Hadamard regularization 
scales. Thus, all the Hadamard regularization scales, one for each particles (they were 
previously denoted s±, • • • , sn in [HI), are to be replaced by the unique scale £q. For instance, 
in the equivalent representation of Hadamard's partie finie as an analytic continuation, 
making explicit the presence of those arbitrary constant scales, we have 

P f(yoH i FP / d 3 X f-r..^r^M ,3,5) 



47T 



where the symbol FP is understood as the finite part of the Laurent expansion of the integral 
when all of the a a 's tend to zero. The Hadamard partie finie ( 13. 24ft or ( 13. 25ft is extremely 
convenient to implement in practical computations. 



2. IR divergence 

Next we have also to worry about the IR-type divergence of the Poisson potential P, due 
to the behavior of the source F at spatial infinity, when r — > +oo. Indeed, we find that the 
near-zone potential X we have to evaluate (and which is to be computed at 1PN relative 
order) is given by an iterated Poisson integral which is divergent at infinity. The appearance 
of a divergent near- zone potential X(y a ) is a novel feature of the present calculation; indeed 
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the problem did not arise in the previous computation of the 3PN equations of motion 
because we needed instead the gradient (diX)(y a ), which is convergent. 

Fortunately, the problem of IR divergences has been solved in the general case, for any 



isolated PN source and at any PN order. Here we shall follow the formalism of Ref. [42 
which uses systematically a regularized version of the Poisson integral which is appropriate 
for solving the hierarchy of PN equations. The idea is to introduce inside the Poisson integral 
a regularization factor (r/r Q ) B , where B is a complex number, and where r is an arbitrary 
IR scale (as £q is an arbitrary UV scale). The regularized Poisson integral is then defined as 
the finite part (FP) coefficient, i.e. the coefficient of the zeroth power of B, in the Laurent 



expansion of the integral when B — > 0. It was proved in [42| that the latter regularized 
Poisson integral is a solution of the Poisson equation for a general regular (smooth) source, 
and is amenable to iteration up to any PN order. 

In principle the latter procedure is defined in 3 dimensions. However since we are here 
solving the <i-dimensional field equations, we shall first define it in d dimensions, so that the 
far-zone part of the Poisson potential reads 

where FP refers to the finite part when B — > 0. The precise meaning of considering the 
FP process on a <i-dimensional integral has been discussed in |Bj]. 15 Here, since we consider 
only the far-zone part of the integral free of UV divergences, we can immediately take the 
limiting case e = and get 

pSV) = --Lfp/" d'x^T^M (3.27) 



4vr J r>n \r 

The result will depend both on the IR cut-off scale tq and intermediate radius 1Z, but we 
shall check that these constants disappear in the final results. At the point y a we have 

Pf (yoH -± F p/ ,^) B £!!1M. (3 . 28 ) 

4vr J r> Ti \r J r a 

Finally, the sum of the near- zone integral (I3.25P and far-zone one (13.281) gives our complete 
prescription for the finite part of the dimensionally regularized Poisson integral as 



We do not detail how this integral is computed in practice but refer to previous works [28 
and, for the treatment of the bound at infinity, Sec. IV C of J43J. We checked that the 
sum of the resulting UV and IR-regularized Poisson integral is independent of the arbitrary 
constant length-scale 71. 



15 It could be possible to use dimensional regularization to cure not only the UV divergences but also the IR 
ones (without the FP when B — > 0). However this would imply major modifications of the PN iteration 
scheme for a general source; this has not been attempted. 
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IV. POST-NEWTONIAN RESULTS 



A. The regularized 3PN metric 



The post-Newtonian metric is generated by a system of two point particles, and computed 
at the location of the particle 1 following the prescription (jl.6p . Here we shall somewhat 
abusively simply denote this metric g a p{yi,t), or in short g a p(yi), and similarly for other 
quantities evaluated at the location of particle 1, so that 



g a p(yi 



AC 



lim^(x,i) 

.x->yi K 



(4.1) 



We compute all the required near-zone potentials V, Vi, • ■ ■ , X, f at point 1 (actually, only 
X and T are new in the present computation with respect to Ref. 29]) and regularize them 
according to the procedure of the previous section. The regularized metric in harmonic 
coordinates is now obtained in closed analytic form up to 3PN order as 16 
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The masses mi and ni2 have an arbitrary mass ratio. As usual we denote by ri2 = |yi — y2 1 the relative 
distance between the two particles in harmonic coordinates, by 1112 = (yi — y"2)A*i2 the direction from 
particle 2 to particle 1, and by V12 = vi — V2 the relative velocity, where v a = dy a /di is the coordinate 
velocity of particle a. The Euclidean scalar product between two vectors A and B is (AB). Parentheses 
around indices are used to indicate symmetrization, i.e. A' l B^' = ^{A l B^ + A^B 1 ). 
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We indicate explicitly the post-Newtonian remainders 0(c n ). This metric agrees up to 2PN 



order with the already known result obtained in 25J, and recently used in 24( for obtaining 
the self-force at 2PN order. Because of the helical Killing symmetry we did not include here 
the 2.5PN radiation-reaction terms; these can be found in Eqs. (7.6) of 25 . 

In some logarithmic terms at 3PN order £ denotes the arbitrary constant length scale 
associated with dimensional regularization, which relates the d-dimensional gravitational 
constant G^ to the usual Newton constant G through (13.41) . This scale appears conjointly 
with the numerical combination 

p=V4^e c/2 , (4.3) 

where G = 0.5772 • ■ ■ is the Euler-Mascheroni constant. 17 

Notice the important feature that the metric in harmonic coordinates involves some poles 
oc 1/e at the 3PN order in the 00 and Oi components, where e is related to the spatial 
dimension d by d = 3 + e (see Sec. IIIip . and formally tends to zero. The results presented 
in Eqs. (14.21) include the pole part ~ e' 1 and the complete finite part ~ e°, and neglect the 
terms tending to zero when e — Y 0; for simplicity we do not indicate the remainders 0(e). 

However there is an exception to the above rule, in that we have to re-introduce the 
correction terms 0(e) in the Newtonian part of the metric. Indeed, when we shall reduce 
the metric to the center-of-mass frame and then to circular orbits, these corrections will 
be multiplied by poles at 3PN order, and will contribute in fine to the finite part at 3PN 
order. Such corrections will be necessary only in the 00 component of the metric, where the 3- 
dimensional Newtonian potential at the location of the particle 1, namely Vn(2/i) = Gm 2 /r 12 , 
is to be replaced by its d- dimensional version 18 
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7 The number p appears in the expansion when e — > of the parameter k defined by Eq. ()3.14j) as k = 
l-e\np + 0{e 2 ). 

8 Terms C(e 2 ) are neglected. See (|A17[) and (|A18[) for the exact expressions of the Newtonian potential 
and acceleration in d dimensions. 
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The poles oc e~ x in the metric (14. 2 j) could be removed by a coordinate transformation and 
a suitable shift of the two particle's world-lines. This is discussed in Appendix lAl where we 
compute the regularized 3PN metric using the alternative Hadamard regularization. There 
we show that, modulo some assumptions necessary to overcome the known drawbacks of 
Hadamard's regularization (viz the presence of ambiguities at 3PN order), the regularized 
metrics in the two regularization schemes are physically equivalent, in the sense that they 
differ by a coordinate transformation plus the additional effect of some shifts of the world- 
lines of the particles. In particular we find complete agreement with the shifts necessary to 
link together the 3PN equations of motion computed in Hadamard 28 1 and dimensional j29[ 
regularizations. However, for the present purpose it is better to leave as they are the poles 
oc e~ l in the metric (14. 2p . because we are going to compute a gauge invariant quantity, and 
the poles will ultimately be automatically cancelled in the final result. 

Finally we note that the metric depends also on the extra arbitrary constant r , present 
in some logarithmic terms of the 00 component of the metric at 3PN order. This constant 
comes from the IR regularization of the metric at spatial infinity, as discussed in Sec. IIII B 21 
and it shall also disappear in the final gauge invariant result. 

As an important check of the metric (14. 2 p we have verified that it is invariant under 
a general Lorentz boost, considered in a perturbative 3PN sense. The Lorentz invariance 
permits checking most of the 2PN terms and also the dynamical 3PN ones. The only terms 
which are not checked by a 3PN Lorentz boost are the 3PN static ones — those that do not 
depend on velocities. 



B. The gauge invariant quantity u 

To compute the gauge invariant quantity u T (associated with particle 1 for stationary, 
circular orbits), we adopt its coordinate form as given by (II. 3p . namely 

= i^-gM^j , (4.5) 

and plug into it the 3PN regularized metric explicitly obtained in ( 14.21) . To begin with, this 
yields the expression of u l at 3PN order for an arbitrary mass ratio q = mi/m 2 , and for a 
generic non-circular orbit in a general reference frame. 

We then choose the frame of the center of mass (CM), which is consistently defined at 
the 3PN order by the nullity of the 3PN center-of-mass integral of the motion deduced from 
the 3PN equations of motion jijj . We want to express the individual positions y a = y„ M 
and velocities v a = v^ M (with a = 1,2 labelling the particles) relatively to the center of 
mass in terms of the relative position y 12 and relative velocity v 12 . We know how to do this 



at 3PN order in Hadamard regularization [45j, and we know that the particle's trajectories 



in Hadamard regularization differ by a shift of world-lines from those computed with dimen- 
sional regularization 29[. So in order to get y^ M and v^ M in dimensional regularization we 



apply directly the shift of world-lines on the known expressions in Hadamard regularization; 
this is detailed in Appendix [B] 

Having replaced the positions and velocities by their CM expressions ya M [yi2,vi2] and 
v^ M [yi2, V12], the quantity u l becomes a functional of yi 2 and vi 2 which we now reduce to 
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the case of circular orbits. This means that (7112^12) = exactly, 19 and that the relative 
orbital velocity squared v\ 2 takes a specific expression in terms of the relative separation r±2 
or, rather, in terms of the particular dimensionless post-Newtonian parameter defined by 
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where m = mi + rri2 is the total mass of the binary. 20 We find in Appendix 
required relation, valid in dimensional regularization, is 
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where ^ = mxm^jvn 2 is the symmetric mass ratio, related to the asymmetric mass ratio q by 
v — q/(l + q) 2 - From now on we assume that ra\ ^ m 2 to prepare the ground for the small 
mass ratio case mi rri2 in which v = q + 0(q 2 ). Notice the presence of a pole oc 1/e at 
the 3PN order in ( 14 .7p . and recall that £0 is the dimensional regularization scale, and that p 
is defined by (I4.3p . Note also that we have included the O(e) correction in the Newtonian 
approximation of the expression (14. 7ft ; this is crucial because multiplying that Newtonian 
term 0(e) by a quantity having a pole at 3PN will yield a finite part contribution at 3PN 
order. The last step of the calculation consists of replacing 7 by its expansion in powers of 
the convenient alternative dimensionless gauge invariant PN parameter x, directly related 
to the orbital frequency Q = v 12/^12 by 
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To find 7 as a power series in x to 3PN order we invert (14.71) and obtain 
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When finally replacing 7 by x we discover most satisfactorily that all the poles oc \ je 
cancel out in the final expression for w*, as well as the associated constant to (and the 
pure number p). Furthermore, the IR constant r also disappears from the result when 
parameterized by the frequency-related parameter x. No matter what the mass ratio, our 
final result for a 3PN, algebraic relationship between u T (to which u l now evaluates) and x 
(or equivalently Q), is: 
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Consistently with the helical Killing symmetry we neglect radiation-reaction effects. 

Recall that the orbital separation 7*12 is here dchncd in harmonic coordinates, and differs from the 

Schwarzschild coordinate distance r used in the SF calculation of Sec. [TTJ 
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where we denote A = (m 2 — m\)jm = yl — 4z/, so that the test-mass limit of particle 1 
corresponds to v — > 0. The expression (14.101) is a polynomial in x with coefficients depending 
only on the symmetric mass ratio v\ it is therefore clearly gauge invariant. While it has been 
shown in |24| (see also Sec. II Bl above) that u T is gauge invariant at any PN order, in the 
extreme mass ratio limit v <C 1, here we find that it is also gauge invariant for any mass 
ratio up to 3PN order. 21 



V. COMPARISON OF POST-NEWTONIAN AND SELF-FORCE RESULTS 



We now reduce the 3PN expression (I4.10p in the small mass ratio regime q = m\jm2 <C 1. 
We express the result in terms of the non-symmetric PN parameter introduced in (I2.16p . 
which is more suited than x to the small mass ratio limit of particle 1, namely 
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Using x = y(l + q) 2 ^ 3 and v = q/(l + q) 2 we obtain, up to say the quadratic order in q 
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This is to be compared with the result of SF calculations, which take the general form 

u T = «L + q «sf + Q 2 "psf + C(q 3 ) , (5.3) 

with self-force and post-self-force coefficients u§ F and Wp SF respectively. From f)5.2p we thus 
recover the 3PN expansion of the Schwarzschildean result, i.e. 
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Next, we obtain the self-force contribution u^ F up to 3PN order as 
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(5.5) 



21 As a test of the initial expression of it* for a generic orbit in a general frame (i.e. before going to the CM 
frame), we checked that dw*/cU = after reduction to circular orbits, as required by the helical symmetry, 
i.e. neglecting the radiation-reaction. 
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FIG. 2: The self-force contribution Ug F to u T plotted as a function of the gauge invariant variable y^ 1 . 
Note that y" 1 is equal to Rn/ m 2, an invariant measure of the orbital radius, scaled by the black hole mass 
7712 [cf. Eq. (|2.16|) ]. The "exact" numerical points are taken from Ref. [24 1. 



The 2PN result is in agreement with (I2.17P as it should. For the much more difficult 3PN 
coefficient, whose value depends on subtle issues regarding the self-field regularization (see 
Sec. IIIII) . we thus find 

C 3 pn = -^ + |vt 2 . (5.6) 

We get also the 3PN expansion of the post-self-force, which could be compared with future 
SF analyses with second-order black hole perturbations, 22 

u T PSF = V + 3y 2 + fy 3 +(J§- f/) y 4 + Otf) , (5.7) 

as well as all higher post-self-force effects up to 3PN order. 

Numerically, the 3PN coefficient in the self-force is C3PN = —27.6879 • • • . This shows a 
remarkable agreement between the post-Newtonian prediction and the result of the numerical 
SF calculation reported in (I2.18p . namely C|p N = —27.677 ± 0.005. The two results are 
consistent at the 2a level with five significant digits. This agreement can also be visualized 
in Fig. [2j where we show the SF contribution «g F to u T as a function of as well as 
the successive Newtonian, 1PN, 2PN and 3PN approximations to u§ F . Observe notably 



Notice that u|p < and Up SF > (at least up to 3PN order) . The effect of the self- force is to reduce the 
value of u T , while the post-self- force tends to increase it. 
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FIG. 3: Numerically derived residuals, i.e., after removal of the 2PN and 3PN self- force contributions to 
w|p, plotted as a function of the gauge invariant variable y^ 1 . Compare with scales in Fig.O Note that y^ 1 
is equal to Rn/rri2, an invariant measure of the orbital radius, scaled by the black hole mass m-i [cf. (|2.16[) ]. 



the nice convergence of the successive PN approximations toward the exact SF result. The 
3PN approximation is roughly 1% accurate up to y~ x = 10, and roughly 5% accurate up 
to y~ l = 7, not very far from the highly relativistic Schwarzschild innermost stable circular 
orbit (ISCO) for which y^ co = 6. 23 

This successful comparison between SF and PN calculations confirms the soundness of 
both approximations in describing compact binary systems. In the post-Newtonian calcula- 
tion, this encompasses the post-Newtonian expansion as applied to the binary equations of 



motion (28|,|47J, and includes the treatment of the issues associated with the UV divergencies 



using dimensional regularization 29|, |30] . In addition, the IR divergences, too, in the PN 
calculation (see Sec. IIHB 2\i are seen to be correctly treated, since their effects vanish in 
the final result (14. 101) . In the perturbative self-force calculation embodied in fl2.17p - fl2.18p . 
this includes the delicate handling of gauge and the numerically taxing split of the metric 
near the particle into singular and regular pieces following the prescriptions in [l3j]. In this 
light, it would be interesting to address the opposite question, namely that of estimating the 
accuracy of the black hole perturbation formalism by comparing several truncated self-force 
series to the "exact" PN result in the slow motion limit. This would require at least a 



2.", 



4a ] for a recent calculation of the shift of the Schwarzschild ISCO induced by the conservative part 



Sec 



of the self-force. 
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second-order perturbative SF calculation. 

Our post-Newtonian calculation contains additional results which have not been used in 
this paper. For example, it already contains some of these higher-order self-force terms, as 
evidenced by (15.71) . Similarly, our numerical self-force calculation actually contains much 
more information than is indicated by the numerical coefficient we give in (j2.18p . This 
is most simply illustrated in Fig. [31 where we show, over the large- range used for our 
numerical fitting, the full 2PN and 3PN residuals, that is, the residuals after the known 
2PN and (now) 3PN terms have been subtracted from our numerical data. In fact, we have 
gone to considerable lengths to ensure that we would have high quality numerical data to 
work with here. The smooth curve of 3PN residuals, several orders of magnitude below 
the 2PN curve, is a testament to this data quality and represents the starting point for an 
investigation which more adequately explains the appropriate higher order PN nature of our 
numerical data; especially the presence of logarithmic terms in higher PN approximations. 
The pressing need for this explanation is strong motivation for further work 37 . 
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Appendix A: Relation to Hadamard regularization 

As an important check of the DR calculation of u T , we have also performed the complete 
calculation using the alternative Hadamard regularization (HR), in the variant proposed in 



28J and called the "extended Hadamard regularization". The HR is essentially based on 
the Hadamard partie finie integral (I3.24p . Unfortunately we know that the HR, even in the 
extended variant, is not entirely satisfying and cannot give a completely consistent picture at 



the 3PN order |28l. 1471] . In particular it yields some ambiguities in the form of a small number 
of arbitrary parameters which cannot be computed within this regularization. However, the 
HR gives the correct answer provided that the ambiguity parameters are specified by some 
external arguments, or fixed by comparison with some non-ambiguous calculations. Then 
it becomes a non trivial check to show that it is possible to adjust a few HR ambiguity 
parameters so that the complete result, which is generally made of many more terms, agrees 
with the result of DR. 

The ambiguity parameters in HR come from the unknown relations between two sets of 
arbitrary length scales denoted s a and r' a (where a labels the particles). Here the scales 



s a are introduced into the Hadamard partie finie [40|, |41[ of Poisson integrals with singular 
sources, when computed at any field point x' different from the singularities y a . The s a 's 
appear when the Poisson integrals develop logarithmic divergences on the singular points 
(i.e. when the source point x over which one integrates equals y a ). The other scales r' a come 
from the singular limit of the Poisson potential when the field point x' itself tends toward 



28 



the singularity y a ; hence we have in fact r' a = |x' — y„| which shows up in the form of some 
"constant" hir' a which is formally infinite. It was shown that the relation between the scales 



s a and r' a must involve the masses m a (and m = X]& m f>)' and is necessarily of the type [28 



r' \ „m 



In ^ = a + (3 — . (Al) 
\SaJ m a 

Here a and (3 denote some purely numerical constants and are called ambiguity parameters. 
After imposing the link flAip to get rid of the scales s a , it was shown that the remaining 
scales r' a are gauge constants which can be removed by a change of gauge. 24 In the case of 
the 3PN equations of motion (EOM), and for the extended variant of HR, it was found that 
the correct values are 29| 25 



159 , „ 1987 
«eom = — and /?eom = -— • (A2) 

In the case of the 3PN mass quadrupole (MQ) moment needed to compute the 3PN radiation 
field, the values using the same extended HR turned out to be [H, 43 



2(> 



9451 

«mq = -9240 and = • ( A3 ) 

The fact that the MQ values are different from the EOM values already shows that the 
HR is not satisfying: Two different computations give inconsistent determinations of the 
ambiguity parameters. For the present computation of the quantity u T and comparison 
with SF calculations, we have shown that the extended variant of HR reproduces exactly 
the result of DR [i.e. f)4.10p above] if and only if we have the still different values 

asF = "l54 and fl^iio" (A4) 

Although this result shows again that the Hadamard regularization is not consistent at 
3PN order, we argue that it constitutes a powerful check of our calculation, because for the 
complete agreement we have to adjust no more than two unknown coefficients. In particular 
we find that the HR constants r' a , which remain after imposing the relation flAip . cancel out 
properly in the final result. The HR calculation is also interesting because it corresponds to 
a different harmonic coordinate system and a different definition of the particle's world-lines. 
Since we have the regularized 3PN metric in both HR and DR, we can now study in detail 
their difference — adopting the values fjA4[) in the HR scheme. 

We shall find that the two metrics differ by an infinitesimal 3PN coordinate transfor- 
mation in the "bulk", i.e. outside the particle's world-lines, and also by an intrinsic shift 
of these world-lines. In particular we shall recover the total shift obtained at the level of 
the equations of motion in Ref. j29[, but shall prove that this shift is made of the latter 
intrinsic shift, plus the shift induced by the coordinate transformation. Let the coordinate 



24 Hence the fact that lnr' a is actually "infinite" docs not really matter. 

25 More precisely, q;eom was determined by requiring that the equations of motion should derive from a 
Lagrangian formulation, but /3eom (which was denoted A in Refs. 28 . 
fixed by DR. 

6 The coefficients «mq and Pmq were respectively denoted £ and k in Refs. [4|,|43(. 



had to wait until its value was 
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transformation between the two metrics be 5x a = e a (x), where e° = 0(c~ 7 ) and e l = 0(c~ 6 ) 
as appropriate to 3PN order. The transformation of the bulk metric is (for simplicity we 
omit the PN remainders) 

SeQoo = -2<9 e - e*0i0oo - _ e \y2)^rr » ( A5a ) 

dy\ dy\ 

S e g 0i = -2d {0 ei) , (A5b) 

6 e gij = -29(i€j) . (A5c) 

The terms involving partial derivatives with respect to the source points y a take into account 
the shifts of the trajectories e(y a ) through a modification of the source dependence of the 
metric g a /3(x;y a ,v a ); we use the notation y a = (ct,y a ) and v a = (c, v a ). Since e a is of order 
3PN, the (?oo's in the RHS of ( lA5aj) are simply Newtonian. At the point 1 we get 27 



d d 
Se9oo(Vi) = -29 e (yi) - e *(yi)-Q-[ [9oo(yi)] - el (y^Q^i [9oo(Vi)] , ( A 6a) 

Se9oi(yi) = -29 ( oei)(yi) , (A6b) 
8e9ij(yi) = -2%e j) (?/ 1 ) . (A6c) 

Now we have found that in order to relate the two metrics one must additionally perform a 
shift n a of the particle's trajectories at the 3PN order, i.e. n a = 0(c~ 6 ). Such shift will be 
"intrinsic" in the sense that it will not be induced by any coordinate transformation of the 
bulk metric. It yields the additional change of the metric components evaluated at point y\\ 

dy\ dy\ 

while there is no change in the other components at that order, i.e. S K goi(yi) = S K gij(yi) = 0. 
Our final result is that the two regularized metrics at point 1 are related through 

9af(yi) = 9af(yi) + SegafiiVl) + &K9c#(yi) ■ (A8) 

This relation is a functional equality relating the two metric functionals at point 1, whose 
coordinates yf can be seen as dummy variables. 

An important check of our finding (|A8[) is that it contains the previous result derived at 
the level of the 3PN equations of motion (not considering the bulk metric), namely that the 
DR and HR equations of motion merely differ by some shifts £ a of the particle's world-lines. 
This result established the physical equivalence of HR and DR at the level of the equations 
of motion {29|. Indeed, we discover that the total shift found here, which is made up of the 
shift e(y a ) induced by the coordinate transformation plus the intrinsic shift K a , is precisely 



equal to the shift of the world-lines of the particles found in [29|. Thus 



L = e(y a ) + n a ■ (A9) 

So we have proved that £ a is partly intrinsic and partly induced by a coordinate transforma- 
tion of the bulk metric. Since there is no physics involved in a coordinate transformation, it 



27 Note that our too compact notation <?oo(2/i) stands in fact for <?oo(2/i; Ua,v a )', thus we have used the obvious 
relations d[g 00 (yi)}/dy\ = {d l g 00 )(yi) + {dg o/dy\)(yi) and d[g 00 (y 1 )]/dyi = {dgm/dy^yx). 
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can be argued that only the intrinsic part of the shift n a corresponds to the physical process 



of renormalization which was performed in [29 



We give now the explicit expressions. The coordinate transformation at any field point 
x = (ct, x) reads 



e (x) = - 



edx) 
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Remind that the DR-related quantites p, £o and k are defined by Eqs. ( 14.31) . ( 13. 4ft and ( 13. 14[) 
respectively, while r 2 is an HR constant of particle 2. The symbol 1 -H- 2 means adding 
the previous expression [i.e. excluding the term /i(t) in (jAlOaft ]. but with all particle labels 
exchanged (thus r' 2 would be changed to r[). Note that r[ and r' 2 are considered as true 
constants in (IA10I) . The function f±(t) is given by 



, u . 91 G 3 m\m 2 k(l + e) . . 
= 77 7 2TT~ yn 12 v 12) 
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Notice that because of the presence of fi(t) in (jATUa]) , the time component of the gauge 
vector 6q(x) is not symmetric by exchange 1 -H- 2. This coordinate transformation satisfies 
the harmonic gauge condition De a = in d dimensions at the 3PN accuracy. We note also 
that the time component eo of the coordinate transformation has a pole part oc s" 1 , as well 
as a pole-independent part, but that the space component e« is pole-free. 

Beware that strictly speaking ( 1A10|) is not the coordinate transformation between the HR 
metric g^{x) and the DR metric g^{x) in the bulk. It is solely the restriction e a (yi) of this 
gauge transformation at the location of particle 1 that correctly relates the two regularized 
metrics g^ivi) an d Q^fijji) at that location. Indeed, if the gauge transformation (lAlOj) 
was to be valid in the bulk, it would induce poles in the 00 and Oi components of g^{x). 



But it was shown in (29|, based on diagrammatic arguments, that such poles, at 3PN order, 



can only be present in the 00 component gffi{x) of the DR metric. However, the restriction 
e a (yi) at the location of particle 1 of the pole-free gauge transformation in the bulk does 
generate poles in the Oi component of the DR metric at y\ [see Eq. (14.2bj) ]. 

The shift induced by this coordinate transformation is pole-free, and we immediately get 
from (lAlObj) [up to a correction 0(e)] 



(A12) 



where is the <i-dimensional Newtonian acceleration of body a given by (1A18|) below. 
Because e(y a ) does not contain any pole, we observe from ( lA6cj) that the spatial part of the 
regularized metric will be free of poles at 3PN order as well [cf. Eq. (I4.2cj) ] . Next we find 
that the additional shift n a does contain a pole, and explicitly reads 
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so that the total shift as defined by (1A9|) is given by 

11 GW 



2 In 
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1540 



(A14) 



in perfect agreement with the result of [29 



For completeness we now give the result for the difference between the two regularized 
metrics, 5g a p{yi) = g^f(yi) - QafiVi)- Combining (TA6])-flAE]) with flA10ll - flA13l) we get 
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The end result for the Hadamard regularized 3PN metric, g^iui), then follows from com- 
bining the previous difference with the explicit expression (14. 2 p for the DR metric. One 
can check while performing the sum that all poles oc 1/e and the associated ^-dependent 
logarithmic terms cancel out, so that the HR result only depends on the UV gauge constants 
r' a and also, of course, on the IR regularization constant r . 

Note that the DR metric (1A8|) is really the metric experienced by the particle in d = 3 + e 
dimensions. It is thus very important to include in that metric all corrections of order e 
which could yield finite contributions after multiplication by quantities involving poles. As 
already mentioned, for the problem of computing u T for circular orbits we have to write the 
Newtonian part of the 00 component of the metric as = — 1 + 2V^ d) /c 2 + C(c- 4 ), where 
the Newtonian potential satisfying the d- dimensional Poisson equation AV^ = — -iirG^a^ 
with Newtonian source density <tn = 2< ^~^ m a 5^ is given by 



2(d - 2) , G^m, 



d 



+?*£ 
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(A16) 



with DR value at point a [see also (14. 4p ] 
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In the same vein the Newtonian acceleration in (1A13|) - (1A14|) should read 



a: = Vtf' W = -^X;^n„, (A18) 



ab 



Appendix B: Circular orbits in d dimensions 

In this Appendix we describe our way to reduce a general rf-dimensional expression such 
as the regularized metric (14. 2ft — valid for arbitrary binary orbits and in a general frame (in 
harmonic coordinates) — to the center-of-mass (CM) frame and then to circular orbits. The 
relevant formulas to do so have been worked out at 3PN order within HR (see jsj for more 
details), and we need here the corresponding formulas valid in DR. Basically we shall rely 
on the HR results and apply to them the known shifts of the particle's world-lines to deduce 
the corresponding DR results. 

The 3PN equations of motion of compact binaries using HR turned out to depend on 



one, and only one, ambiguity parameter called A (denoted /?eom in Appendix fA} [281] . a nd 



2S 



to be physically equivalent to the DR equations of motion if and only if A = [29| 
This means that the difference between the DR and HR accelerations of body 1 (say) is 
exclusively due to a shift of the world-lines of the particles y a — > y a + £ a through 



.DR HR 



af H | A __!987 + ^ai . (Bl) 

1 3080 S 



The explicit A-dependent expression of the 3PN-accurate acceleration can be found in 
Eq. (7.16) of j28j. The effect of the shifts £ a on the acceleration of body 1 is 

8 a N 

*ai = - & + 0(c~ 8 ) , (B2) 

where £J 2 = £! — £2, and the dot stands for a derivative with respect to coordinate time t. 
The shift £ a has been given in (1A14|) above; recall the presence therein of a pole oc e~ x . To 



be consistent one needs to include in the Newtonian acceleration a^ the corrections of order 
e, and the correct expression to do so is given by ( 1A18I) . 



By definition, the CM frame is such that the center of mass position G vanishes. Within 
HR, we have G HR = when the individual positions of the particles y a are given as some 
functionals of the relative position y 12 and velocity v 12 according to y a = y™[yi2,v 12 ]. 
The explicit expression of the functionals ya R [yi2, V12] up to 3PN order can be found in 
Eqs. (3.6)-(3.7) of [iol ]. Similarly, within DR we shall have G DR = when the individual 
positions of the particles are related to the relative position and velocity according to some 
new functional relations 

Va = y° R [yi 2 ,v 12 ], (B3) 

which we want to determine. Now, by the effect of the shifts of the world-lines the expression 

of the center of mass position in DR will be different from that in HR, and be given by 
g dr = qhr + where 

<5 ? G = -m^i - m 2 & + 0(c~ 8 ) . (B4) 



28 This result is equivalent to the one of Ref. 30|; see also [48|,|49| for an alternative, ambiguity- free derivation 
of the 3PN equations of motion. 
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Therefore we find that the DR functionals ya R [yi2,vi2] are related to the HR functionals 
y^ R [yi2, V12] through y° R = y™ + S^y a 29 with the same shift for both particles given by 



5& a 



nil 



m 



m 2 



m 



6 + o{c 



(B5) 



The DR expressions (1B3|) are thus easily determined from the HR results. 

Next, from the DR equations of motion flBll) -f )B2|) in a general frame, we go to the CM 
frame by replacing the individual positions and velocities by the relative ones according to 



y a = Ya R [yi2, V12] and also v a 



y? R [yi2,vi 2 j 



Turning off the well-known 2.5PN radiation- 



reaction terms, and restricting the result to circular orbits [thus (2/12^12) = 0], we get the 
relative acceleration of the binary within DR in the form af 2 R = — Q 2 yi2, where the orbital 
frequency Q can then be computed iteratively as an expansion in powers of the PN parameter 
Gm/(ri 2 c 2 ), with rn = \~yn\- To 3PN order we find 
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As in the Hadamard case, we invert this relation to express 7 as a PN series in powers of 
the gauge invariant parameter x = (Gmf2/c 3 ) 2 / 3 , with result 
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(B7) 



x 3 + 0{c 



The equations (1B6I) - (1B7I) are the DR equivalent of Eqs. (188) and (191) in [8|, which are 
valid in the coordinate system used in the HR case. Of course the results coincide up to 
2PN order as they should. Notice however that we kept the terms proportional to e in 
the Newtonian terms of (1B6|) - (IB7|) . because in the process of computing u T the Newtonian 
terms will get multiplied by some poles e~ x occurring at 3PN order, and these corrections 
will contribute to the final result. We are done for the results necessary for the computation 
of u T for circular orbits as reported in Sec. IIVI 

As a useful check, we compute the total energy of the binary for circular orbits within 
DR, making use of Eqs. flB5j) -( lB7|) . For a generic orbit and in a general frame, the DR 
energy functional is related to the HR one through E BR = E RR + S^E. The effect of the 
shifts £ a on the energy explicitly reads 



5i-E = -mi v\ Ci - m 2 v\ Q + f j 2 



where the Newtonian gravitational potential energy in d dimensions is 

{d) _ 2(d-2) Ggmimg 
d-1 r? 2 " 2 ' 



(B8) 



(B9) 



Note that we mean by this a functional equality, valid for any dummy variables yi2 and via- 
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At this stage, we use the expression of the total energy E nR as computed within HR, and 
given e.g. by Eq. (170) of [8], and add to it the term 5^E defined by flB8l) - (1B9j) . Our first 
check is that this 3PN-accurate energy E DR for a generic orbit in a general frame within 
DR is conserved, i.e. E DR = when neglecting the 2.5PN radiation-reaction terms. This 
requires consistently order reducing the result, i.e. replacing the accelerations in the time 
derivative of E DR using the DR equations of motion (lBip - (IB2l) . 

Now, we obtain the expression E DR in the center-of-mass frame by replacing the individual 
positions and velocities by their expressions ya R [yi2,vi2] and y° R [yi2, V12]. Restricting 
ourselves to circular orbits, the resulting CM energy depends only on vf 2 = ?"i2^ 2 an d 7. 
Then we replace vf 2 by its PN expansion in powers of 7 using (IB6I) . and finally replace 7 by 
its PN expansion (IB7I) in powers of x. We find that all poles oc e^ 1 disappear in the process; 
therefore we can take the limit e — > 0, and get the gauge invariant expression 

^ nR mvc 2 f / 3 v\ / 27 19 v 2 \ 2 

E = x < 1 + )x + + — v )x 2 

2 1 V 4 12/ V 8 8 24/ 
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675 
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which coincides with the well-known 3PN expression of the total energy for circular orbits 
as given e.g. by Eq. (192) of 0. 
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